Enhanced and reduced atom number fluctuations in a BEC splitter 
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We measure atom number statistics after splitting a gas of ultracold *^Rb atoms in a purely 
magnetic double-well potential created on an atom chip. Well below the critical temperature for 
Bose-Einstein condensation Tc, we observe reduced fluctuations down to — 4.9dB below the atom 
shot noise level. Fluctuations rise to more than +3.8 dB close to Tc, before reaching the shot noise 
level for higher temperatures. We use two-mode and classical field simulations to model these results. 
This allows us to confirm that the super-shot noise fluctuations directly originate from quantum 
statistics. 
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Since the achievement of Bose-Einstein condensation 
(BEC) in dilute atomic gases, different experimental 
techniques have been developed in order to coherently 
split a BEC into two spatially separate parts [IHl], 
with atom interferometry as one of the motivations. 
Even though BECs are usually in the weakly interacting 
regime, the interactions between the particles dramati- 
cally affect the physics of the splitting. In particular, 
repulsive interactions limit the phase coherence between 
the two split parts [F , but also reduce atom number dif- 
ference fluctuations, giving rise to non-classical squeezed 
states [SHin]- 

In this Letter, we use a purely static magnetic poten- 
tial created on an atom chip to realize a nonlinear spatial 
"beam" splitter for a BEC. We investigate the physics of 
the splitting and focus on atom number fluctuations and 
the role of temperature. At low temperatures, where the 
interaction energy dominates, we directly observe num- 
ber squeezed states with relative population fluctuations 
—4.9 dB below shot noise, as first shown in [TU: and indi- 
rectly observed in |S]. The two separated but weakly 
linked parts of the BEC constitute a bosonic Joseph- 
son junction, usually described by a two mode model 
(TMM) 11 . Our results are in agreement with the 
TMM, which also predicts that the observed squeezing 
is accompanied by high phase coherence. 

The magnetic trap configuration allows barrier heights 
up to several /iK and straightforward evaporative cool- 
ing, so that we can separate clouds with increasing tem- 
perature all the way to the non-degenerate regime. In 
the intermediate temperature regime, where both a sig- 
nificant condensate and thermal fraction are present, we 
observe large super-binomial fluctuations in the number 
difference between the two parts. This excess of fluctua- 
tions is a direct signature of the Bose statistics, in close 
analogy to the bunching effect in quantum optics [16j . 



Close to the BEC transition, the condensates show sig- 
nificant depletion and the TMM breaks down. We com- 
plement our experiments by a theoretical investigation of 
this regime using a classical field approach and show that 
large super-binomial fluctuations are a general feature at 
thermal equilibrium. Although the experiments are not 
performed at equilibrium, our observations are still in 
qualitative agreement with these theoretical results. 

Our experiment uses a two-layer atom chip to pre- 
pare a *^Rb BEC in the |F = 2, m-p = 2) hypcrfinc state, 
and then split it in a double-well potential with ad- 
justable barrier height and well spacing. The magnetic 
potential along the weak axis is well approximated by 
V{x) = Vb[l — (x/xq)'^]'^, where 14 is the barrier height 
and 2x0 is the distance between the wells (Fig. [l] (a) 
and (b)). The transverse potential is harmonic with trap 
frequencies tUy^z/^Tr w IkHz. Unlike RF dressed poten- 
tials, which tend to create very elongated traps ^ iS]j 
our traps have an aspect ratio close to 1, which strongly 
suppresses phase fiuctuations within each well. We split 
the trap by increasing I2 , which simultaneously increases 
both Vb and xq (Fig. [I]^b)-(e)). The interesting regime 
of two weakly coupled condensates occurs for a barrier 
height on the order of the chemical potential (few kHz), 
corresponding to a current I2 — 2.5 mA in our experi- 
ments. In this region, the tunneling rate is divided by 2 
when the barrier Vb/h increases by 75 Hz. We estimate 
fiuctuations of the barrier height to be on the order of 
20 Hz. A magnetic field component normal to the chip 
leads to a small energy difference AE between the two 
potential minima. A magnetic shield with an attenuation 
factor of ^ 30 reduces fiuctuations of AE to a few Hz. 

The expected atom number fiuctuations are on the or- 
der of a few tens of atoms, placing stringent requirements 
on the imaging system. We perform absorption imaging 
using a back-illuminated CCD camera and 3 ms delay be- 



(a) d = 4|jm 



yJ^Lji 



L 



(b) 


li 


2 ll 


10 


"1 


1 


8 


-1 


i 


W 6 


-1 


/ 




\ 


/ 


~4 


_ \ 


/ 


CQ 


\ 


/ 


2 


- \ u. 


i / 


n 


L iViX V 


N/, 1 



-4 4 




I, (mA) 



FIG. 1: (color online), (a) Chip geometry: we use DC 
currents superposed with a static homogeneous field Bo = 
{Ba:, By, Bz) = (15, 8, 0) G. To first approximation, the poten- 
tial can be understood as follows: The field of Jo = 100 mA 
combines with Bo to create harmonic confinement in the xy 
plane. The minimum is located at a distance 20 = fioIo/^TrBy 
from the 7o wire and the field in the minimum is B^. This 
minimum is modified by 7i = 2.0mA and I2, whose fields are 
predominantly along x on the trap axis. For zq > d, the two li 
wires together provide harmonic confinement along x. 72, of 
opposite polarity, creates the barrier of adjustable height and 
also determines the spacing between the two resulting wells, 
(b) Profile of the trapping potential along the splitting axis 
X, for I2 = mA (cooling trap, dashed line) and I2 = 2.4 mA 
(solid line). The x axis is the same as in (a), (c) Barrier height 
Vb (top), trap frequencies {fx, fy, fz) (center) and position xo 
of the right minimum (bottom) as functions of the current I2 . 
Lines: numerical calculation based on current configuration. 
Blue circles: measured frequencies fx- 



tween absorption and reference images. The total quan- 
tum efficiency including optical losses is q — 0.84. Probe 
pulses have r = 50 /iS duration and an intensity close 
to saturation. Atom numbers are calibrated following a 
procedure inspired by [12]. We estimate our total sys- 
tematic error to be at most 18%. Atoms in the two wells 
are resolved by applying a large gradient with the central 
conductor (72) for the first 50 /xs of a 6 ms time-of-flight. 
Photon shot noise, scaling as the root of the image area, 
leads to a standard deviation of about 20 atoms for a 
BEC of 10'^ atoms imaged over 800 pixels, and about 
60 atoms for a thermal cloud (2 x 10^ atoms, 7000 pix- 
els), pixel area being 9.5 /im^ in the object plane. This 
is close to the theoretical limit ^yl6/qa^T = 0.19 /im^^, 
where a is the scattering cross section and F the natu- 
ral linewidth, and always remains below atom shot noise; 
furthermore, this noise is well-characterized and stable so 
that it can be subtracted. Technical fluctuations (moving 
fringes) constitute a second noise source, which is negli- 
gible for small BECs, but becomes comparable to photon 



shot noise for large atomic clouds. We do not attempt 
to remove this noise, but indicate an estimate of its level 
(green dots in Fig. [s]) . 

To measure the fluctuations of the atom number dif- 
ference iVi — N-[i between the left and the right well, 
we acquire a dataset consisting of a large number of 
absorption images (typically > 100) taken under iden- 
tical conditions. In each image, we measure Nj^ and 
Nr (see inset of Fig. [s]). To compensate for a slight 
imbalance in the splitting, we determine the probability 
to be in the left (right) well Jl^r = {Nl,r/N), where 
N = Nl + Njf and the average is taken over all images. 
For each image, we calculate n = JrNl — IlNr, which 
is the deviation of {Nl — Nb)/2 from the expected value 
(/i ^ fR){NL + Nii)/2. We deflne the number squeez- 
ing factor as the variance of z = nj ^ ji^j^ij^i^ -I- iV/j). 
Correcting for the photon shot noise contribution Zp [TS] 
leads to the final expression ^ = {z^) — (z^). This defi- 
nition is first-order insensitive to fluctuations in the total 
number of atoms and produces ^^ = 1 for a binomial 
distribution. In our data, the correction (Zp) is always 
smaller than (z^) itself. 

In a first experiment, we split an almost pure BEC of 
1300 atoms and investigate the influence of the time Tr 
during which the barrier is raised by increasing I2 from 
zero to 3.9 mA, well above the chemical potential. The 
final number squeezing ^^ is shown in Fig. ^ We observe 
a decrease of fluctuations below the shot noise limit with 
increasing r^ up to 50 ms, followed by a slow increase 
for longer times. This agrees with the expectation that 
an optimum should exist between very short t^ creat- 
ing excitations in the BEC and/or not leaving sufficient 
time for tunnelling, and very long r^ where heating and 
atom loss become important. For a quantitative analy- 
sis, we start by noting that low-energy excitations from 
the symmetric many-particle ground state correspond to 
populating the first excited mode, which is the longitu- 
dinal dipole mode for Vb = and becomes the Josephson 
plasmon mode as the barrier raises. These excitations are 
spaced by hujp, where the plasma frequency Wp decreases 
from LOp — ujx to as the barrier raises [S] . The interest- 
ing dynamics can be expected to occur when the barrier 
height approaches the chemical potential. In this region, 
hujp is already much smaller than hiOx, which is the en- 
ergy scale of the next higher modes. Therefore, a TMM 
is expected to work well in this region at least as long 
as ksT < hhJx- The two parameters entering the TMM 
are the charging energy Ec accounting for the interac- 
tion between the particles and the Josephson energy Ej 
that characterizes the tunnelling. We calculate them for 
each barrier height by solving the 3D Gross-Piteavskii 
equation for the experimental trap [11]. In the TMM 
simulations, we describe the initial state before splitting 
by a thermal density matrix, which we evolve according 
to the von Neumann equation. We start the simulation 
at I2 = 1.9 mA, slightly below the splitting point. We 
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FIG. 2: (color online). First experiment: Fluctuations of 
the atom number diflference as a function of splitting time r^ 
for a BEG with N=1300. Red squares: measured ^^; dotted 
line: dynamical TMM simulation. Upper inset: Ec (dashed 
line) and Ej (solid line) from numerical 3D GPE solution for 
1300 atoms using the full calculated potential. Lower inset: 
Degradation of the number squeezing under the influence of 
heating. The BEG is held for a variable time r^ before split- 
ting; Tr = 50 ms is constant. 



have checked that the results depend only weakly on the 
starting point (initial I2) in this region. An initial tem- 
perature Ti = 50 nK reproduces the measured fluctua- 
tions after a ramp of 10 ms. The Ti thus found is then 
used for all longer t^. The simulation reproduces well the 
observed squeezing for r^. up to 30 ms. For r^ > 40 ms, 
the experimental data shows a degradation of squeezing 
which we attribute to technical heating. To confirm this 
hypothesis, we hold the BEC for a variable time r^ before 
splitting it with a 50 ms ramp (inset of Fig. [2]). 

In a second experiment, we investigate the effect of 
the temperature on the squeezing. We ramp I2 linearly 
from to 3.9 mA in 50 ms, and vary the temperature 
by changing the final frequency i^f of the evaporative 
cooling ramp. The results are plotted in Fig. [3J Below 
10.73 MHz, i.e. well below Tc, we observe a crossover 
from sub- to super-binomial fluctuations with increas- 
ing temperature. Around 10.73 MHz, the condensed 
and thermal fractions are comparable, and we observe 
large super-poissonian fluctuations with a maximum of 
-1-3.8 dB. As the temperature increases further, fluctua- 
tions decrease, and level when the condensate fraction 
reaches zero. The measured asymptotic level of 1 dB is 
consistent with binomial statistics given our fringe noise 
(Fig. IsFa) , green circles) . According to the TMM calcula- 
tion, our best result of ^^ = —4.9 J;q;4 dB corresponds to 
a phase coherence (cosLp) ~ 0.93 immediately after the 
splitting, where if is the relative phase between the two 
clouds. This would result in a possible metrology gain of 
—4.4 dB compared to the standard quantum limit using 
this state in an atom interferometer. 
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FIG. 3: (color online). Second experiment: Effect of the 
temperature on ^^. We vary the final frequency Uf of the 
evaporative cooling ramp before splitting the cloud using a 
50 ms linear current ramp. Typical absorbtion images are 
shown in the inset of (a). Red squares: 5^. Green circles: 
fringe noise. (Estimation of the ^^ we expect to measure 
when the actual fluctuations are binomial - the estimate is 
obtained by analyzing a region of the image that does not 
contain atoms.) (b) Total atom number A^ (red diamonds) 
and condensed fraction Nq/N prior to splitting (blue circles). 



The observed sub- and super-binomial regimes origi- 
nate from the interplay between interactions and quan- 
tum statistics. Lowering the temperature, the onset of 
super-binomial fluctuations occurs when quantum de- 
generacy becomes important. Fluctuations are given by 
the probability distribution of the macroscopic configu- 
rations with a given atom number difference n. This 
distribution is binomial in the classical gas regime, lead- 
ing to (n^) = N/A. In the degenerate regime, the en- 
tropy effect which favors small number differences van- 
ishes and, if each Fock state with a given n is equiprob- 
able, (n^) is as high as N'^/12. The crossover from 
super- to sub-binomial fluctuations comes from the inter- 
action energy cost associated with number fluctuations, 
which eventually exceeds the available thermal energy. 
This low temperature regime is well described by the 



TMM, which predicts fluctuations 



= kBT/Ec, as 



a result of the equipartition theorem. We heuristically 
extend the TMM to higher temperatures by adding a 
binomial contribution for the thermal cloud, leading to 
^2 ^ AkBT/{NEc) + 1 - N2m/N, where N2m is the con- 
densate fraction, i.e. the population of the two lowest 
modes, which has to be determined independently. This 
formula implies the existence of a maximum: on one hand 
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FIG. 4: (color online). Classical field simulations (open sym- 
bols+dotted lines) and modified TMM (full symbols+lines). 
Blue lines and triangles: A'' = 6000; red lines and circles: 
A^ = 17000. The trapping potential is harmonic transver- 
sally {fy,z ~ 1120, 1473 Hz), and quartic along x with Vb/h = 
2.69 kHz and xo ~ 3.8 /im. The longitudinal oscillation fre- 
quency within each well is /j; — 411 Hz. The temperature 
unit in the figure is the ideal gas transition temperature in 
the harmonic potential prior to splitting T^^'^'^ (JV = 6000: 
yharm ^ Q.SQ ^K, N = 17000: T^'^"'" = 1.25 mK). Inset: Frac- 
tion of the population in the two lowest energy modes N2m , 
extracted from the classical field simulations. 



the decrease of temperature increases A^2mj which con- 
tributes to the super-binomial signal. On the other hand, 
interactions that tend to lower the fluctuations dominate 
more and more as the temperature decreases. Given 
that the maximum occurs at some fraction of the crit- 
ical temperature T^., where a macroscopic population in 
the two lowest mode appears, this maximum scales as 
ksTc/iNEc), which is large for a cold gas in the weakly 
interacting regime. 

In order to account more accurately for the thermal 
cloud contribution, we perform a multimode calcula- 
tion in the classical field approximation. The fluctua- 
tions are decomposed into a shot noise term N plus a 
term where the fields appear in the normal order. We 
compute the latter by sampling the Glauber-P distribu- 
tion that we approximate by the classical distribution 
P ex exp{-f3E[i(j,il;*]}, where E[iI;,iIj*] is the Gross- 
Pitaevskii energy functional [TS]. We include the first 
quantum correction to the classical field, which has the 
effect to change the shot noise term N into TV — N2m- As 
expected, the classical field approximation predicts an 
increase of fluctuations starting above T^, when the gas 
enters the degenerate regime. The sub-binomial regime 
is out of reach of the classical field approximation us- 
ing the Glauber-P distribution. The results are shown in 
Fig. |4J and qualitatively reproduce the experimental re- 
sults (which were obtained in a dynamical process) . The 
inset shows N2m extracted from the classical field calcu- 
lations, which we also use in the TMM curve. In both 



models, fluctuations are largely independent of the atom 
number N for a given condensed fraction, justifying the 
comparison with the experimental data. 

These results show the interplay of quantum statis- 
tics and interactions in a simple and fundamental finite- 
temperature system. The double well can be seen as 
a "nutshell" version of the Bose-Hubbard model whose 
more complex dynamics also lead to the Mott insulator 
transition [17] . Our results also highlight the generic fea- 
tures of spatial splitters for trapped BEGs: interactions 
allow for the creation of non-classical states of poten- 
tial interest for quantum metrology. However, the same 
interactions will lead to phase spreading with a rate pro- 
portional to ^vN Ec after the splitting. Hence, tuning 
the interactions to reduce Ec after splitting is a neces- 
sary step for the use of such a non-linear splitter in a real 
interferometer. 
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